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ABSTRACT 

The  design  and  implementation  of  hardware  for  approximate 
digital  division  is  described.  Using  standard,  commercially- 
available  T  L  components,  a  compact  divider  has  been  devel¬ 
oped  which  is  well-suited  for  use  within  5  MHz  systems  hav¬ 
ing  three-state  busses.  After  forming  double-length  denom¬ 
inator  reciprocals  by  table  look-up  to  high  accuracies, 
single  and  extended-precision  quotients  are  made  available 
by  consecutive  numerator  multiplications  using  a  single¬ 
chip  array  multiplier. 
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INTRODUCTION 


I  . 


Tor  a  large  number  of  digital  signal  processing  algo¬ 
rithms,  it  has  generally  been  true  that  the  infrequent  need 
for  a  full-precision  2's  complement  divide  capability  has  re¬ 
resulted  either  in  this  function  being  implemented  in  software, 
or  else  denominator  reciprocal  values  being  prestored  in  memory 
for  later  add-shift  numerator  multiplication  to  achieve 
"real  time"  division.  For  those  applications  wherein  a  fast, 
compact  hardware  divide  capability  is  essential,  the  choices 
have  been  somewhat  limited.  For  example,  one  may  implement 

either  a  successive-approximation  or  a  non-restoring  division 

2 

algorithm  using  one  or  two  handfuls  of  MSI  T  L  devices  or 
utilize  a  single-chip  multiplier/divider  (e.g.,  the  MMI  #67516), 
hut  in  either  case  the  expected  quotient  calculation  time  exceeds 
two  microseconds. 

A  table  look-up  method  is  described  in  this  paper 

for  forming  the  approximate  reciprocal  of  any  16-bit  2's 

2 

complement  input  using  12  MSI  T  L  chips  (plus  SSI  gates) .  Ac¬ 
cess  time  is  slightly  more  than  100  ns.  When  used  together 
with  a  single-chip  16  x  16  array  multiplier,  single-precision 
quotient  calculations  require  about  250  ns. 

Figure  1  illustrates  how  the  scaled  reciprocal  function 
fK/X)  may  be  used  within  a  two-bus  qC  structure  (K  =  2  *5). 


INST  No.  ACTION 

1  K  /  X  * —  X  ,  M  U  LT  *• —  Y 

2  MULT  —  K/X,  Y  —  MULT 

3  MULT  —  K/X,  Y  —  MULT 


Fig.  1.  Two-bus  microcomputer  structure  showing  data  flow 
for  single  and  double -precision  ' pseudodivides ' . 


In  the  context  of  a  5  MIPs  microcomputer  currently 
being  developed,  nominal  200  ns  instruction  times  provide  generous 
overhead  allowances  for  data  bus  transfers  and  control  delays, 
while  expanding  the  total  alloted  time  for  reciprocal  estimation 
and  multiplication  to  400  ns.  All  functional  blocks  shown  in 
Figure  1  communicate  via  three-state  data  buses  X  and  Y,  where 
K/X  is  equipped  with  an  input  register,  the  multiplier  is  fully 
registered  on  both  input  and  output  ports,  and  the  data  source/ 
sink  is  intended  to  represent  one  of  several  possible  two-ported 
data  elements,  such  as  a  data  memory  or  CPE  register  file.  After 
operands  X  and  Y  are  transferred  and  K/X  is  formed  during  in¬ 
struction  1,  a  single-precision  quotient  is  generated  during 
instruction  2  by  loading  K/X  into  the  multiplier,  doing  the 
multiplication,  and  delivering  the  result  to  the  data  sink. 
Instruction  3  accesses  the  extended-precision  portion  of  the  K/X 
table  for  multiplication  by  the  same  numerator  Y  to  complete  the 
formation  of  an  approximate  double-precision  quotient. 

When  the  reciprocal  table  is  used  for  the  "pseudodivide" 
operations  described  here,  scaling  factor  K  is  eliminated  by 
exploiting  the  output  select  feature  (MSP/LSP)  of  the  multiplier 
to  rescale  the  resultant  output,  Q  =  Y/X.  A  transfer  function 
such  as  K/X  whose  hardware  is  internally  registered  and  equipped 
with  a  three-state  port  may  be  applied  within  single-bus  structures 
as  well  as  dedicated  (non-programmable)  hardware. 
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II.  K/X  APPROXIMATION 

There  is  ample  incentive  to  minimize  the  required  amount 
of  table  storage  space  for  a  reciprocal  function  which  will  ac¬ 
commodate  16-bit  signed  2’s  complement  inputs.  Even  a  single- 

20 

precision  (16-bit)  unsampled  table  requires  2  bits  of  storage 
or  sixty- four  16K-bit  programmable  read-only  memory  devices, 
apart  from  support  circuitry.  By  imposing  a  reasonable  error 
constraint,  we  can  piecewise-approximate  the  K/X  curve  and  get 
a  dramatic  table  memory  (PROM)  saving.  Thus,  assuming  for  now 
that  output  data  precision  is  not  an  issue  (i.e.,  assuming 
"infinite  precision") ,  we  seek  to  approximate  K/X  with  closely 
spaced  points  in  the  region  of  steepest  slope  (smallest  X) ,  but 
to  use  progressively  more  sparse  spacing  of  points  as  the  slope 
diminishes  and  ]X|  approaches  one.  For  the  moment,  let  us  con¬ 
sider  only  the  positive  region  of  the  function  R  =  K/X.  Figure  2 
shows  the  alternate  point  sampling  that  would  be  likely  within 
a  region  of  steeper  slope. 

Reduction  of  the  required  number  of  table  entries  in¬ 
volves  1)  the  choice  of  compromise  K/X  values,  and  2)  the  choice 
of  progressively  wider  K/X  sample  spacing  as  X  increases  from 
decimal  zero  to  one.  A  PROM  table  addressing  algorithm  is 
needed  which  will  compress  an  unsampled,  2^-word  table  into 
a  much  smaller  one.  For  logical  simplicity,  let  the  address 
quantization  be  binary-weighted. 


Fig.  2.  Illustration  of  alternate  point  sampling 
for  R  =  K/X  CK  =  2-15,  x  =  n-2'15). 


A 


For  the  positive  region,  assuming  2's  complement  frac¬ 
tional  inputs,  define  index  n  =  0-*-32767  ,  X  =  n-2 and  K  =  2 
such  that 

32767 

R  =  K/X  =  215/n- 2~15  =  1/n  n=Q  (1) 

If  the  quantization  factor  (A)  represents  the  number  of  con¬ 
tiguous  X  points  spanned  for  a  given  table  address,  denotes 

/\ 

normalized  approximation  error,  and  R  is  the  approximate  value 
for  R  over  the  same  span,  then  for  uniform  normalized  peak 
error , 

Ea  =  (1/n-R) /  (1/n)  =  [R- 1/  (n+A  - 1)]  /  [l/  (n+A  - 1;]  (2) 

from  which 

R  =  2/[ 2n+ (A- 1)]  (3) 

and 

Ea  =  (A- 1) /[2n+ (A- 1)]  (4) 

Reciprocal  function  K/X  can  be  organized  into  several 
neighboring  zones,  each  having  its  own  unique  A.  As  X  =  n‘2  * 
increases  from  n  =  1,  this  first  zone  (#1)  requires  that  A  =  1 
(from  (4),  A  =  2  would  imply  that  EA  =  0.33  at  n  =  1,  an  unac¬ 
ceptably  large  error) .  Zone  1  is  somewhat  special  in  that 
A  =  1  results  in  a  1:1  relationship  between  table  addresses 
and  data,  and  no  approximation  error  exists.  The  end  point 
for  zone  1  may  be  deduced  by  calculating  the  beginning  value 
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of  n  for  zone  2,  subject  to  an  error  constraint.  From  (4),  for 
Ea  <  (A-l)/[2n+(A-l)]  ,  n  >  (1 -EA)  (A- 1)/2EA  (5) 

Thus  if  we  require  EA  _<  1/21,  then  n  >_  100.  Using  the  nearest 
larger  binary  value,  choose  n  =  128  as  the  starting  point  of 


zone  2, 


To  minimize  the  number  of  table  entries  as  X  =  n-2 


increases  toward  1,  we  wish  to  increase  A  to  2A  at  the  earliest 
opportunity.  For  uniform  peak  error,  vhe  zone  boundary  can  be 
determined  by  solving  (4)  for  n  when  A  =  2A,  where  n^  and 
n^+i  are  the  first  address  points  in  zones  i  and  i+1,  respec¬ 
tively  : 

Ea  =  (A.-l)/[2ni+(A.-l)]  =  (2A.-l)/[2ni+1  + 

(  2  A  4  -  1  j]  (6) 

(ni  +  P/ni  =  (2Ai-l)/(Ai-l)  (7) 

The  starting  point  of  zone  3  (a  =  4)  determines  the  end  point 
of  zone  2,  and  from  (7)  is  n^_^  =  128  (4-l)/(2-l)  =  384. 

Equation  7  determines  the  beginning  points  of  subsequent  zones 
in  a  similar  manner,  with  uniform  peak  error  throughout,  pre¬ 
determined  by  the  starting  point  of  zone  2.  Because  zone  2 
begins  at  n  =  128,  universal  EA  =  0.00389.  Table  1  shows 
the  range  of  n  values  for  the  9  zones  needed  to  generate  the 
positive  portion  of  an  approximate  K/X  function.  Fortunately, 

A 

zone  9  may  be  merged  with  zone  8  by  letting  the  last  stored  R 
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TABLE  1 

NUMBER  RANCH  AND  QUANTIZATION  FACTORS  (A) 
FOR  A  9-ZONE  APPROXIMATE  RECIPROCAL  FUNCTION 


value  span  n  =  32512  -*•  32767  ,  A  =  256.  Each  zone  requires  128 
storage  locations,  and  the  8- zone  arrangement  reduces  the  needed 
storage  by  a  factor  of  32  with  respect  to  unsampled  storage. 

For  both  positive  and  negative  regions  of  the  function  K/X,  with 
approximation  error  of  less  than  1 / 2 %  and  a  data  wordlength  of 
16  bits,  32K  bits  of  PROM  table  storage  is  sufficient. 

Figure  3  shows  the  function  R(x)  (at  somewhat  distorted 

/v 

scale)  for  which  the  approximate  function  R(x)  is  to  be  com¬ 
puted  and  stored  in  PROM  tables.  The  intent  here  is  to  illus¬ 
trate  the  numerical  range  of  interest,  and  to  take  note  of  de¬ 
partures  from  the  ideal  curve  attributable  to  2's  complement 
asymmetry  along  with  a  limit  when  X  =  0.  For  the  positive 
region,  X  ranges  from  0  to  1-2  with  R(0)  and  R(2 

constrained  to  the  largest  representable  positive  fraction, 

- 1 5  - 1 5 

1-2  .  For  the  negative  region,  X  ranges  from  -2  to  -1, 

and  R(x)  is  the  straightforward  scaled  reciprocal.  Excluding 
R(0)  and  R(2  *  ),  the  asymmetry  about  the  origin  may  be  ex¬ 
pressed  as 

R(A+)  =  - R (A+- 1 )  »  (8) 

where  A+  is  any  given  address  within  the  positive  region. 
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1  I  [ .  K/X  QUANTIZATION 

It  is  important  to  consider  the  accuracy  limits  of 
K  =  K/X,  i.e.,  the  precision  with  which  one  may  represent  the 
scaled  reciprocal.  Let  us  define  the  normalized  quantization 
error  as 


Eq  =  (R-Rn)/R, 


Q  ‘ 


f9) 


where  R  is  the  exact  (infinite  precision)  value  of  K/X  and  Rg 
is  its  quantized  value.  As  indicated  in  Fig.  4,  the  absolute 
amount  of  quantization  error  in  a  single-precision  R  table  may 
approach  the  value  of  the  LSB  (2  1^).  But  on  a  percentage 
basis,  this  error  has  the  least  impact  when  the  absolute  value 
of  X  is  small,  i.e.,  when  R  is  large. 

Tables  2  and  3  indicate  the  expected  worst-case  total 
errors  (approximation  plus  quantization)  for  single  and  double¬ 
precision  realizations,  respectively.  Recalling  from  (1)  that 

I 

max 

(10) 


R  =  K/X  =  1/1 


max 
1  =  0 


where 


I  •  2 


-15 


(ID 


both  tables  show  normalized  error  as  a  function  of  I  and 

max 

X  together  with  the  input  values  (I  ,  and  X  ,)  at  which 
max  p  k  p  k 

total  error  ET  is  maximum.  As  expected,  Table  2  shows  approxi- 
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Fig.  4.  Maximum  quantization  error  (Eq)  versus  denominator 
(x)  magnitude  for  single-precision  reciprocal,  R. 
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mation  error  (E^)  to  be  zero  for  zone  1,  and  constraining  X  to 

-  8 

be  no  larger  than  2  results  in  quantization  error  (Eq)  of 
less  than  0.4%  as  indicated  in  Fig.  4.  Note  also  that  while 
approximation  error  is  minimal  throughout  the  usable  input 
range,  quantization  error  rapidly  becomes  unacceptable  with 
increasing  X.  However,  for  those  applications  requiring  good 
accuracy  over  the  full  dynamic  range  of  X,  Table  3  shows  that 
extended  precision  approximation  of  K/X  yields  improved  results. 
Simply  stated,  a  s ingle -precis ion  approximation  to  K/X  can  re¬ 
solve  to  within  2  ^  (the  LSB)  ,  whereas  a  double-precision 
approach  is  accurate  to  with  2  The  penalty  for  extended 

precision  is  an  additional  multiply  operation  when  doing  a 
"pseudodivide" . 

The  Fortran  programs  for  generating  Tables  2  and  3  are 
shown  in  Appendix  A. 
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For  single-precision  results  to  have  tolerable  limits 
on  error,  some  constraints  must  be  imposed  on  the  range  of 
allowable  divider  inputs  shown  in  Fig.  4.  if,  after  compensat 
ing  for  scaling  we  say  that  16-bit  quotient 

Q  =  Y/X  ,  (12) 

and  avoidance  of  overflow  requires  that 

| Y |  <  ] X |  ,  (13) 

then  for  quantization  error  Eq <  0.41  as  chosen  for  this  par¬ 
ticular  design, 


00 

I 

Cs] 

V| 

X 

(14) 

Q 1  >  2'7 

(15) 

In  the  context  of  a  given  pC  program  requiring  a  fast 
divide,  a  priori  knowledge  and  maintenance  of  the 
relative  magnitudes  of  X  and  Y  can  permit  straightforward, 
fast  generation  of  single-precision  quotients.  If,  however, 
a  number  of  instructions  are  expended  in  first  shifting 
divider  operands,  the  hardware  divide  feature  may  have  less 
appeal.  For  this  reason,  the  K/X  address  map  function  has 
been  designed  to  accommodate  the  full  denominator  (X)  dynamic 
range,  and  permits  the  addressing  of  a  double-precision 
(32-bit)  K/X  table  simply  by  including  two  additional  PROMs. 

/N 

The  divider  hardware  then  provides  a  single-precision  Q  in 

A 

two  instructions  (400  ns),  and  an  extended  precision  Q  with 


one  additional  200  ns  instruction,  subject  only  to  the  initial 
requirement  for  |Yj  <  |X|.  If,  on  the  other  hand,  we  had  elected 
to  limit  the  input  range  to  |Xj  _<  2  (see  Fig.  4),  a  simple 
256  x  16  PROM  table  with  no  address  mapping  would  have  sufficed. 
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IV. 


ADDRESS  MAPPING 


Figure  5  is  a  simplified  block  diagram  showing  the 
generation  of  K/X  by  cascading  address  map  and  table  memory 
sections.  The  principal  design  task  here  is  the  derivation  of 
a  mapping  function  whose  implementation  will  have  reasonable 
speed  and  a  minimal  number  of  integrated  circuits.  By  examining 
the  2's  complement  representation  of  X  as  described  in  Table  1, 
the  logical  relationships  for  address  mapping  may  be  deduced 
and  minimized.  Although  only  numbers  for  the  positive  region 
are  considered  here,  the  same  mapping  logic  may  be  used  for  the 
negative  region  after  inverting  input  word  X  under  control  of 
sign  bit  X^.  Figure  5  shows  an  11-bit  table  address,  where 
bit  A^q  determines  the  table  region  (+,-),  address  bits  A^ 
through  Ay  control  zone  selection,  and  bits  A^  through  AQ 

constitute  the  zone  address. 

Let  seven  intermediate  variables  be  defined  as  follows: 


C 

C 

C 

C 

C 


1 

2 

3 

4 

5 


C 


6 


C 


7 


x8x7  *  x8x7 
x9x8x7  *  x9a8  *  X7) 

X10X9X8X7  +  X10^X9  +  ^8  +  X7  ^ 

X11X10X9X8X7  +  X11^X10  +  X9  +  ^8  +  X7  ^ 
X12X11X10X9X8X7  +  X1 2  ^ X1 1  +  X10  +  X9  +  X8  +  X7^ 


X13X12Y11X10X9X8X7  +  X13^X12  +  X11  +  X10 

x9  +  Xg  +  X?) 


X14X13X12X11X10X9X8X7 


(16) 
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ADDRESS 

MAP 


TABLE 

MEMORY 


Simplified  reciprocal  function  generator 


Then 


A10  Xs 


(]  7) 


and  zone  select  bits  are 


A9  X14X13X12C4  +  X14X13CS  +  X14C6 

+  x14x13x12x11x10xgx8x7  +  x14 


A8  =  X14X13X12X11X10C2  +  X14X13X12X11C3 
+  X14C6  +  X14X13X12X11X10X9X8X7  +  X14 


A7  "  X14X13X12X11X10X9C1  +  X14X13X12X11C3 
+  X14X13C5  +  X14X13X12X11X10X9X8X7  +  X14 


(18) 


The  required  zone  address  bits  are  shown  in  Table  4  below. 


TABLE  4 

ZONE  ADDRESS  BITS  AS  A  FUNCTION  OF  ZONE 
NUMBER  AND  INPUT  VARIABLES  X,C 
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The  foregoing  expressions  are  the  result  of  examining 


I 

I 


the  relationships  between  binary  input  values  (X)  and  required 
address  values  (A),  with  particular  emphasis  upon  the  X  states 
which  mark  the  zone  boundaries.  Simply  stated,  knowledge  of 
the  zone  within  which  any  input  X  falls  determines  1)  thL-  page 
in  table  memory  containing  the  appropriate  K/X  value,  and  2)  the 
modulus  by  which  changes  in  X  advance  the  address  count.  Con¬ 
sider,  for  example,  the  bottom  row  of  table  4.  As  we  proceed 
along  the  K/X  curve  from  zone  1  to  zone  8,  the  table  address 
(with  LSB  =  Aq)  is  incremented  modulo  1,2, 4, 8,...  etc.,  as 
expected . 


I 


I 

i 


! 


V. 


IMPLEMENTATION 


The  schematic  diagram  of  Fig.  6  is  the  result  of  consid¬ 
ering  several  approaches  to  minimization  of  the  address  mapping 
logic,  using  standard  components.  While  intuition  initially- 
suggested  the  use  of  barrel  shifters  or  FPLAs ,  the  use  of  fast 
PROMs  in  conjunction  with  programmable  multiplexers  yielded 
the  most  compact  design.  The  zone  address  pattern  of  Table  4 
is  suggestive  of  a  large  multiplexer,  and  indeed  the  PMUXs  of 
Fig.  6  permit  generation  of. all  zone  address  bits  except  A^. 
PROMs  A  and  B  provide  all  the  variables  of  Eqs .  16  and  18,  plus 

A^ .  The  16K-bit  table  storage  PROMs  represent  the  best  density/ 
speed  combination  currently  available,  and  provide  needed 
three-state  outputs.  The  25LS2519  input  registers  are  particu¬ 
larly  well  suited  to  this  application  because  of  their  output 
invert  feature.  While  the  address  mapping  logic  is  configured 
to  accept  binary  inputs  in  magnitude  format,  its  usefulness  is 
extended  by  the  input  register's  ability  to  supply  the  l's 
complement  of  negative  inputs  under  sign  bit  control.  This 
same  sign  bit  selects  the  table  region  (+,-)  via  address 
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One  additional  hardware  trick  is  needed  to  realize  the 
zone  8,9  column  of  Table  4.  As  mentioned  earlier,  a  residual 
zone  9  (32640  £  n  £  32767)  exists  for  which  intermediate  vari¬ 
able  Cy  is  uniquely  true.  This  small  region  can  be  merged  with 
zone  8  with  no  discernible  error  penalty.  In  Table  4,  the 
presence  of  the  Cy  terms  within  the  zone  8,9  column  is  the 
assertion  of  zone  9.  However,  because  PMUX  input  ports  are 
already  fully  utilized  (so  much  so  that  is  generated  via 
PROM  B  instead  of  PMUX  A  or  B) ,  another  method  for  including 
the  contribution  of  Cy  must  be  used.  This  is  done  by  storing 
zone  select  bits  (Ag  ->  Ay)  within  PROM  A  in  their  complemented 
form,  and  applying  Cy  to  PROM  A  as  a  chip  enable.  The  program¬ 
ming  of  PMUXs  and  table  PROMs  is  then  done  in  accordance  with 
Ag,S'g,S’y  instead  of  Ag,Ag,Ay,  and  whenever  Cy  is  asserted,  zone 
selection  is  constrained  to  zone  8  via  the  PROM  A  chip  enable. 
Zone  address  bits  Ag  through  A^  are  unaffected. 

A.  PROM  and  PMUX  Programming 

Referring  to  Fig.  6,  PROM  A  is  required  for  zone  selec¬ 
tion  within  the  table  PROMs,  whereas  PROM  B  is  used  principally 
for  the  development  of  intermediate  address  mapping  variables. 
With  the  exception  of  bit  A6 ,  PMUXs  A  and  B  generate  the  table 
address  shown  in  Table  4.  The  requisite  data  patterns  to  be 
stored  were  written  in  accordance  with  Tables  1  and  4  and 
Eqs.  16  and  18,  and  are  shown  in  Appendix  B.  For  purposes  of 
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electrical  programming,  the  296 9 3  PMUX  is  viewed  as  a 
256  x  4  PROM. 

From  Fig.  6,  Table  PROMs  4  and  3  are  selected  for  single¬ 
precision  readouts,  while  PROMS  2  and  1  contain  extended- 
precision  data.  A  Fortran  program  called  TABLF1  (Appendix  C) 
calculates  the  decimal,  integer,  and  hex  values  required  for 
the  single-precision  PROM  pair,  based  on  parameters  such  as 
number  of  zones  and  the  quantization  factor  within  each  zone. 
TABLE!  outputs  four  disk  files  entitled  BYTE43  DEC,  BYTE43 
INT,  BYTE 4  HEX,  and  BYTE3  HEX.  The  first  two  outputs  pro¬ 
vide  fractional  and  integer  decimal  representations  of  the 
entire  single-precision  list  for  checking  purposes.  Because 
these  files  are  somewhat  voluminous,  they  are  not  included  in 
the  appendix.  The  remaining  two  files  are  the  result  of 
splitting  the  integer  file  into  two  pieces  and  adopting  the 
hexadecimal  format,  such  that  the  resultant  tables  are  directly 
useful  as  data  files  for  programming  the  two  byte-wide  memories, 
PROMs  3  and  4.  These  listings  are  shown  in  Appendix  D. 

Fortran  program  TABLE2  (Appendix  E)  is  nearly  the  same 
as  TABLE1,  except  that  its  purpose  is  to  generate  the  data 
patterns  needed  for  extended-precision  PROMs  2  and  1.  It  does 
so  by  calculating  to  high  precision  all  the  table  values  needed 
for  double-precision,  then  deducting  the  values  already  stored 
in  PROMs  4  and  3,  and  using  the  residue  as  data  for  PROMs  2 
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1 

and  1.  The  TABLE2  output  files  are  labeled  BYTE 2 1  DEC,  | 

BYTE21  INT,  BYTE2  HEX,  and  BYTE1  HEX.  The  latter  two  are  ' 

shown  in  Appendix  F.  As  a  point  of  interest,  program  TABLE2 

resolves  each  table  entry  to  within  2  ^  instead  of  2 

because  a  sign  bit  must  be  included  within  both  the  most  and 

least  significant  16-bit  words  made  available  on  the  K/X  data 

bus . 

The  array  multiplier  to  be  used  in  conjunction  with  the 
reciprocal  function  unit  for  doing  'pseudodivides '  is  the  TRW 
model  MPY-16HJ.  Figure  7  shows  the  block  diagram  for  this 
single-chip  LSI  device  which  can  typically  provide  32-bit 
products  within  125  ns.  Two  particularly  useful  features  in 
the  present  context  are  fully-registered  inputs  and  outputs, 
plus  the  availability  of  three-state  outputs.  Further  refer¬ 
ence  to  these  features  will  be  made  when  defining  multiplier 
operating  modes. 

B.  Reciprocal  Function  Modes 

Control  of  the  reciprocal  function  unit  involves  some 
simple  choices  relating  to  1)  the  status  of  the  three-state 
port  (input  vs  output) ,  2)  whether  or  not  the  data  register 
stores  new  information,  and  3)  single  or  extended-precision 
output  (MSR  or  LSR) .  Table  5  below  outlines  some  basic  modes 
based  upon  the  availability  of  an  input  clock  and  control 
bits  Rq,  R^ ,  and  OE^. 
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CLK  X 


RS 

FT 
CLK  M 


Fig.  7.  TRW  MPY-16  HJ  16  x  16  bit 
parallel  multiplier  block  diagram. 


TABLE  5 

ESSENTIAL  CONTROL  STATES  FOR  RECIPROCAL  FUNCTION  UNIT 


QEr 

R1 

R0 

CLK 

CE43 

CE2 1 

Mode 

1 

0 

0 

0 

1 

0 

Read  MSR 

1 

1 

0 

0 

0 

1 

Read  LSR 

0 

X* 

0 

0 

0 

0 

Hold 

0 

X 

1 

t 

0 

0 

Calc/Hold 

Efl 

don '  t 

care 

As  shown  in  Fig.  6,  Rq  gates  an  external  clock  (C)  and  con¬ 
trols  updating  of  the  data  register.  R.^  determines  the  MSR/ 
LSR  choice  when  accessing  the  table  PROMs  in  conjunction  with 
bus  control  OER.  Referring  to  the  three  instructions  of 
Fig.  1  for  forming  approximate  quotients,  the  corresponding 
reciprocal  unit  commands  in  accordance  with  Table  5  are  Calc/ 
Hold,  Read  MSR,  Read  LSR. 

C.  Multiplier  Modes 

Given  a  set  of  externally- appl ied  clocks  and  control 
states,  the  internal  registers  and  three-state  capability  of 
the  MPY-16HJ  may  be  used  to  full  advantage  for  calculating 
approximate  quotients.  Table  6  defines  several  useful  multi¬ 
plier  modes.  The  use  of  separate  input  and  output  clocks 
(C^  and  C  )  is  assumed,  three-state  ports  become  outputs  when 
TRIM  or  TRIL  are  zero,  and  X  denotes  "don't  care"  conditions. 
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TABLE  6 

CONTROL  STATES 

FOR  TRW 

MULTIPLIER 

CLKx 

CLK 

y 

CLK. 

1  ,m 

TRIM 

TRIL 

Mode 

C- 

l 

0 

C 

0 

1 

1 

Load  X 

0 

c. 

i 

Co 

1 

1 

Load  Y 

c. 

1 

c. 

i 

Co 

1 

1 

Load  X , Y 

X 

X 

0 

1 

1 

Hold 

X 

X 

0 

0 

1 

Read  MSP 

X 

X 

0 

1 

0 

Read  LSP 

Multiplication  occurs  on  any  Load  instruction,  and  the 
resultant  product  is  retained  in  the  output  register(s).  The 
Hold  condition  retains  data  in  the  output  registers,  and  is 
generally  used  when  no  multiplier  outputs  are  driving  the  bus 
The  Read  operations  simply  amount  to  suppressing  the  output 
register  clocks  and  enabling  the  appropriate  three-state 
driver  to  access  either  the  most  or  least  significant  product 

Table  6  depicts  only  a  subset  of  the  total  number  of 
unique  operations  possible  using  the  MPY-16HJ  multiplier.  An 
underlying  assumption  of  Fig.  1  is  that  some  of  the  modes  of 
Table  6  may  be  concatenated  within  any  given  200  ns  nominal 
microcomputer  instruction  cycle.  In  particular,  for  the 
three- instruction  loop  of  Fig.  1,  the  multiplier  does  Load  Y, 
Load  X/Read  LSP,  and  Load  X/Read  LSP  operations  while  the  re¬ 
ciprocal  unit  is  doing  Calc/Hold,  Read  MSR,  and  Read  LSR. 
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The  split  multiplier  instruction  Load  X/Read  LSP  is  faster 
than  others  such  as  Load  Y/Read  LSP  because  the  Y  bus  is  not 
shared  between  input  and  output  traffic,  and  no  three-state 
output  enable  delay  is  incurred.  Another  multiplier  feature 
which  may  be  exploited  is  the  feedthrough  provision.  If  local 
storage  of  a  product  (in  the  multiplier  output  register)  is 
not  needed,  but  instead  a  result  is  to  be  stored  at  some  re¬ 
mote  point,  the  FT  line  on  the  MPY-16HJ  renders  the  output 
register  transparent,  making  the  multiplier  output  available 
earl ier . 

When  doing  division,  the  scale  factor  (K)  associated  with 
the  single-precision  reciprocal  function  is  effectively  removed  by 
shifting  the  multiplier  output  by  2+1^,  i.e.,  by  interpreting 
the  LSP  in  accordance  with  MSP  binary  weights.  For  extended 
precision  calculations,  the  reciprocal  table  output  is  inter¬ 
preted  as  being  unsealed. 

VI.  RESULTS 

Using  commercial  grade  integrated  circuits,  bench  tests 
have  shown  the  clock- to-output  throughput  delay  of  the  one  re¬ 
ciprocal  function  built  to  be  no  greater  than  110  nanoseconds 
at  25°C.  The  expected  delay  value  over  the  commercial  tempera¬ 
ture  range  should  not  exceed  165  ns.  For  those  instances  where 
X  has  been  stored  in  the  K/X  holding  register  for  a  long  time 
(e.g.,  when  fetching  an  extended-precision  result),  access  via 
table  PROM  chip  enable  may  take  as  little  as  50  ns. 
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Tests  over  a  0°  to  70°C  temperature  range  indicate  that 
the  MPY-16HJ  multiplier  typically  requires  125  ns  to  calculate 
32-bit  products.*  The  worst-case  data  sheet  specification  for 
non-pipelined  multiplies  is  175  ns. 

In  light  of  the  above,  the  nominal  200  ns  interval 
cited  earlier  for  forming  reciprocals  and  for  doing  multiplies 
continues  to  be  a  reasonable  estimate. 

Beyond  questions  of  operating  speed,  several  test  cases 
have  been  run  to  verify  the  numerical  validity  of  the  basic 
design,  with  particular  emphasis  on  the  effects  of  approxima¬ 
tion  and  quantization  errors.  Appendix  G  discusses  the  results 
of  two  test  cases,  including  such  issues  as  table  memory  zone 
selection,  scaling,  and  error  calculations. 

The  integrated  circuit  placement  diagram  of  Fig.  8  is 

indicative  of  the  physical  size  of  the  divider.  The  layout 
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shown  requires  about  12  in.  of  circuit  board  area. 

Total  D.C.  power  dissipation  for  the  divider  using 
available  componentry  is  slightly  in  excess  of  8  watts.  This 
number  can  be  reduced  to  5  watts  in  the  future  by  employing 
SPROMs  which  have  a  fast  power-down  capability.  Raytheon 
currently  offers  some  of  the  needed  devices  which  have  a  built- 
in  power  switch,  and  can  potentially  save  70%  of  PROM  power. 

*The  MPY-16HJ  has  been  thoroughly  tested  at  several  temperatures 
using  all  possible  input  operand  combinations  as  part  of  a 
separate  characterization  effort. 
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VI I  . 


SUMMARY 


Tor  those  applications  where  fast,  approximate  division 
is  essential,  the  reciprocal  table  look-up/multiplication 
approach  to  generating  quotients  is  a  viable  alternative.  The 
design  described  in  this  paper  is  of  tractable  size  (15  DIPs) 
and  power  (8  watts) ,  and  provides  double-precision  outputs 
having  total  errors  no  greater  than  0.4%.  If  produced  in  sig¬ 
nificant  quantities,  the  device  count  could  be  further  reduced 
through  the  use  of  semi-custom  integrated  circuits,  particularly 
for  the  address  mapping  function.  Better  accuracies  may  be 
obtained  at  the  expense  of  additional  table  storage.  For  each 
halving  of  the  approximation  error,  twice  as  many  PROM  table 
entries  are  needed. 

The  hardware  implementation  shown  requires  less  than 
200  ns  for  reciprocal  formation  and  for  each  subsequent  multi¬ 
plication.  This  speed  regime  is  well  matched  to  the  capabili¬ 
ties  of  bipolar  (LSTTL)  microcomputers  or  to  dedicated  hard¬ 
ware  systems  running  with  5  MHz  throughput  rates. 
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APPENDIX  A 

FORTRAN  PROGRAMS  FOR  GENERATING  TABLES  3  AND  4 


REAL*8  R  f  RA  f RAQ f  EA  f  ET  f  EQ  f ERROR f EPEAK » XPEAK » EAP r ETP r EQP 
REAL*8  DELfRDELfXI fX2 
WRITE  <  A  *  1 > 

1  FORMAT <//f IXf 'ERROR  ANALYSIS  FOR  APPROXIMATE  1/X'f/> 

WRITE(6f2) 

2  FORMAT (IXf 'DATA  TABULATED  RE  MAX  INPUT  X  -  SINGLE  PRECISION '*/ > 
WRITE ( 6  f  3 ) 

3  FORMAT (IXf  ' EA f EQ f ET=APPROX f QUANT  f  &  TOTAL  WORST-CASE  ERRORS'  f//> 
WRITE ( 6  r  4 ) 

4  FORMAT ( 3X  f  ' IMAX ' f  4X  f ' IPK ' f  6X  f ' XMAX '  f  9X  f 'XPK' f9Xf 'EA' f8Xf 'EQ' ,8X> 
1 'ET' ,//) 

DO  10  L=1 r 16 
IM=2**(L-1 ) 

XM=<  IM--1  ♦  )*2.**<-15) 

EF’EAK --1 . 

11  =  1 

12=128 

J=1 

IMAX=2**(L-1 ) 

IEND=128 

IF  < IMAX . LT . IEND ) IEND=IMAX 
DO  20  1=1 fIEND 
IF(I.GT.2)G0  TO  21 
R=1 .-<2. #*(-15) ) 

GO  TO  22 

21  R=1 ./( 1-1 ) 

22  RA=R 

CALL  QUANT (RAf RAQ) 

GO  TO  50 
20  CONTINUE 

IF( I *EQ. IMAX) GO  TO  40 
15  J=J+1 

13=2* ( 12-11 ) +1 

11=12+1 

12=11+13 

IF ( IMAX  «LT  *12)12=1 MAX 
IF< 12. EQ. 32640) 12=32768 
DO  30  1=11 f 12 
R=1 •/( 1-1 > 

DEL=2.**( J-l ) 

IF ( I .GE«32513)DEL=25A. 

IDEL=DEL 

RDEL=1 ♦ /DEL  ERR0R16  FORTRAN 

X1=I/DEL 

IX1=I/IDEL 
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X2=X1-IX1 

IF  <  X2 » NE .  RDEL ) GO  TO  30 
RA=2./(2.*<I~1 >+DEL~l . ) 

CALL  QUANT (RAfRAQ) 

GO  TO  50 
30  CONTINUE 

IF( I «EQ# I MAX) GO  TO  40 
GO  TO  15 

50  EA=<R~RA)/R 
ET=(R-RAQ)/R 
EQ-ET-EA 
ERROR=ET 

IF  <  ET  *  LT , 0 ) ERROR=-ET 

IF  <  ERROR .  LE  ♦ EPEAK ) GO  TO  51 

EPEAK=ERROR 

IPEAK=I 

XPEAK  =  I*2.**<--15> 

EAP=EA 

ETP=ET 

EQP=EQ 

51  IF <  J . EQ ♦ 1 ) GO  TO  20 
GO  TO  30 

40  UR ITE < 6 » 4 1 ) IM , IPEAK , XM » XPEAK » EAP r EQP > ETP 

41  FORMAT  <lX»I6rlX»I6flXfF12«9»lX*F12*9*lX*F9,6»lX»F9.6*lX>F9,6) 
10  CONTINUE 

STOP 

END 

SUBROUTINE  QUANT < RA , RAQ > 

REAL*8  RA  t RAQ  *  RS » RR  r  C 
DIMENSION  CO ( 16 ) 

RS=RA 
RC^l  ♦ 

IF<RS*G£»0)G0  TO  60 
RC55-!  ♦ 


RS=~RS 

60  C=0 

DO  61  N«2  ? 16 
C0<N)--=2.**<-<N  -1)  ) 
IF(RS.LT*CO(N) )G0  TO  62 
RS=RS-CO(N> 

C*C+CO(N) 

62  RAQ-RC*C 

61  CONTINUE 
RETURN 
END 


ERR0R16  FORTRAN 


t 
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RE AL*8  R » RA , RAQ  *EA*ET*EQ* ERROR * EPEAK * XPEAK » EAP * ETP » EQP 
REAL*8  DEL * RDEL » XI  *  X2 
WRITE(6* 1 ) 

1  FORMAT <//* IX* 'ERROR  ANALYSIS  FOR  APPROXIMATE  l/X'»/> 

WRITE (6*2) 

2  F0RMAT(1X* 'DATA  TABULATED  RE  MAX  INPUT  X  -  DOUBLE  PRECISION'  »/) 
WRITE  (6*3) 

3  FORMAT (IX* 'EA*EQ* ET=APPROX *  QUANT * &  TOTAL  WORST-CASE  ERRORS'*//) 
WRITE (6*4) 

4  FORMAT ( 3X r ' IMAX ' *4X  * ' IPK ' *  6X  * ' XMAX ' *  9X  *  ' XPK' *9X* 'EA'»8X» 'EQ'»8X* 
l'ET'*//) 

DO  10  L=1 *16 
IM=2**(L-1) 

XM=( IM-1 . )*2.**(-15) 

EPEAK=-1 . 

11  =  1 

12=128 

J=1 

IMAX=2**(L-1) 

IEND=128 

IF ( IMAX.LT . IEND) IEND=IMAX 
DO  20  1=1* IEND 
IF< I .GT .2)60  TO  21 
R=l.-(2.**(-15) ) 

GO  TO  22 

21  R=1 . /( 1-1 > 

22  RA~R 

CALL  QUANT(RA*RAQ) 

GO  TO  50 
20  CONTINUE 

IF(  I .EG. I MAX) GO  TO  40 


15  J=J+1 

13=2* ( 12-11 )  +  l 
11=12+1 
12=11+13 

IF ( IMAX « LT . 12 ) 1 2= I MAX 
IF( 12. EQ. 32640) 12=32768 
DO  30  1=11*12 
R=1 ,/( I~1 ) 

DEL=2.**( J-l) 

IF (I .GE.32513)DEL=256. 

IDEL=DEL 

RDEL=1 ./DEL 

X1=I/DEL 

IX 1=1/ I DEL 


ERROR32  FORTRAN 
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X2=X1'-IX1 

I F ( X2 . NE • RDEL >  GO  TO  30 
RA=2  .  / <2 ♦#  < I--1 )  + DEL-1 «  > 

CALL  QUANT <RA»RAQ> 

GO  TO  50 
30  CONTINUE 

IF  < I . EQ ♦ IMAX ) GO  TO  40 
GO  TO  15 

50  EA=(R~RA)/R 
ET=<R~RAQ>/R 
EQ=ET- EA 
ERROR “ET 

IF  <  ET . LT . 0  > ERROR- -ET 
IF(ERROR.LE.EPEAK)GO  TO  51 
EPEAK=ERROR 
I PEAK- I 

XPEAK=I*2.**<-15> 

EAP=EA 

ETP=ET 

EQP=EQ 

51  IF( J»EQ» 1 )G0  TO  20 
GO  TO  30 

40  URITE(6f41)IMf IPEAKfXMf XPEAK f  EAP  f EQP f ETP 

41  FORMAT  <lX»I6rlXflA»lXfF12»9flX»F12«9)»lXrF9*6flXrF9«6»  1XfF9*6 ) 
10  CONTINUE 

STOP 

END 

SUBROUTINE  QUANT (RAfRAQ) 

REAL*8  RAfRAQfRSfRRfC 
DIMENSION  CO (32) 

RS=RA 
RC=1  . 

IF (RS . GE ♦ 0 ) GO  TO  60 
RC—1. 

RS=-RS 

60  C-0 

DO  61  N=2f32 
C0(N)=2.**(-(N~1>) 

IF(RS.LT.CO(N) )G0  TO  62 
RS=RS-CO(N) 

C=C+CO(N) 

62  RAQ=RC*C 

61  CONTINUE 

RETURN  ERROR32  FORTRAN 

END 
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APPENDIX  B 

HEXADECIMAL  DATA  FILES  FOR  PROM  A,  PROM  B, 
AND  PMUX  PROGRAMMING 


PMUX  DATA 

7777777677777765 
7777765377776537 
77765377776537 77 
7653777765377777 
537777773 7777777 
OOOOOOOOOOOOOOOO 
0000000000000000 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 


PROMA  DATA 

7665555444444443 

3333333333333332 

2222222222222222 

1 

1111111111111-111 

1111111111111111 

1111111111111111 

1111111111111110 

OOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOO 
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PROMB  DATA 


I 

I 


00804040COC020202020EOEOEOE01010 
101050505050BOBOBOBOFOFOFOF00808 
08084848484828282828686848689898 
9898D8D8D8D8B8B8B8B8F8F8F8F80404 
04044444444424242424646464641414 
14145454545434343434747474748C8C 
8C8CCCCCCCCCACACACACECECECEC9C9C 
9C9CDCDCDCDCBCBCBCBCFCFCFCFC0202 
02024242424222222222626262621212 
1 2125252525232323232727272720 AO A 
0A0A4A4A4A4A2A2A2A2A6A6A6A6A1 A1 A 
1 A1 A5A5A5A5A3A3A3A3A7A7A7A7A8686 
8686C6C6C6C6A6A6A6A6E6E6E6E69696 
9696D6D6D6D6B6B6B6B6F6F6F6F68E8E 
8E8ECECECECEAEAEAEAEEEEEEEE:E9E9E 
9E9EDEDEDEDEBEBEBEBEFEFEFEFE0000 
0000404040402020202060606060101.0 
10105050505030303030707070700808 
08084848484828282828686868681818 
18185858585838383838787878780404 
04044444444424242424646464641414 
141 4545454543434 3434747474 740C0C 
0C0C4C4C4C4C2C2C2C2C6C6C6C6C1C1C 
1C1C5C5C5C5C303C3C3C7C7C7C7C8282 
8282C2C2C2C2A2A2A2A2E2E2E2E29292 
9292D2D2D2D2B2B2B2B2F2F2F2F28A8A 
8A8ACACACACAAAAAAAAAEAEAEAEA9A9A 
9A9ADADADADABABABABAFAFAFAFA8686 
8686C6C6C6C6A6A6A6A6E6E6E6E69696 
9696D6D6D6D6B6B6B6B6F6F6F6F68E8E 
8E8ECECECECEAEAEAEAEEEEEEEEE9E9E 
9E9EDEDEDEDEBEBEBEBEFEFEFEFE8181 
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APPENDIX  C 

FORTRAN  PROGRAM  FOR  CALCULATING 
SINGLE-PRECISION  TABLE  ENTRIES 


INTEGER*2  YINT < 2048 > 

L0GICAL*1  YPR0M<32> 

EQUIVALENCE  <  YINT  <1 > » YPROM ( 1 ) ) 

REAL* 16  Y < 2048 ) » CF » C  < 16) >  XI » XK » DEL » R 
C(l>=16128. 

C(2)*-192. 

C<3)--4256. 

C(4)=-4240. 

C<5)=-3208. 

C<6)=-2180. 

C<7>=-1410. 

C<8>=~897. 

C  <  9 )  =  -114943 . 

C<10)=~65727. 

C<ll)=-37023. 

C ( 12 ) =-20623 . 

C< 13) =-11399. 

C<  14)  ==-6275. 

C(  15)=s-3457. 

C(  16)~--1920. 

CF=2. 

M-l 

DO  20  J=1 » 16 
IF  <  J  .GE.  9)  M=2 
L2=128*J 
L1=L2 -127 
L3=8*M-J 
DO  20  I -LI *L2 
XI  =  I 

DEL=2.**L3 
XK»C< J)+XI*DEL 

IF( I .EG.897  .OR.  I.EQ.898)G0  TO  15 
IF< I .GT. 897  .AND.  I.LE.1024)G0  TO  10 
IF( I .GE ♦ 1921  .AND.  I.LE.2048)G0  TO  40 
IF( I .EQ. 128  .OR.  I .EG. 1 152>G0  TO  30 
IF (M.EQ. 2 )CF=-2« 

Y(I)=CF/(2.*XK+DEL-1 ,  ) 


TABLE  1  FORTRAN 


GO  ro  20 
10  Y( I )~1 . /XK 
15  Y<897)=1 .  -(2. **(-15) ) 

Y ( 898 )=Y<897) 

GO  TO  20 
40  Y ( I ) a— 1 . /XK 
GO  TO  20 

30  Y<I)=CF/<2.*XK+255. ) 

20  CONTINUE 

DO  70  1-1,2048 
WRITE  <  3, 3  >  Y  < I ) 

3  FORMAT <F36. 33) 

N=1 

ICUM=0 
R=Y ( I ) 

IF( Y< I ) ,GE .0)G0  TO  25 

N=-:l 

R=*-Y  < I ) 

25  R-R-2 •  #*<-- 15 ) 

IF<R.LT.O)GO  TO  35 
ICUM-ICUM+1 
GO  TO  25 

35  YINT<I)=N*ICUM 
70  WRITE <4,4 ) YINT < I ) 

4  FORMAT (18) 

DO  60  1=1,4065,32 

WRITE  (1,2)  YPROM  <  I  > ,  YPROM  ( I  +2  > ,  YPROM  ( 1+4 ) ,  YPROM  ( 1+6 ) , 
lYF'R0M<I+8> ,YPROM <1+10 ) ,YPROM< 1+12) »YPROM< 1+14) , 

1YPR0M(  1  +  16)  ,  YPROM <  1  +  18)  ,YF'ROM(  1+20)  ,YF‘R0M<I+22) , 

1 YPROM ( 1+24) , YPROM ( 1+26) , YPROM ( 1+28) , YPROM <1+30 ) 

60  WRITE < 2, 2 )  YPR0M<I  +  1) ,YPRaM<I  +  3> ,YPR0M<I+5) , YPROM (1+7) » 
1 YPROM < 1+9) , YPROM <1+11 ) , YPROM < 1+13 ) , YPROM < 1+15) , 

1 YPROM< 1+17) , YPROM <1+19 ) ,YPROM< 1+21) ,YPROM< 1+23) , 

1 YPROM < 1+25 ), YPROM < I +27 ), YPROM <1 +29 ), YPROM < I +31 ) 

2  FORMAT <32Z2) 

STOP 

END 


TABLE  1  FORTRAN 
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APPENDIX  D 

HEXADECIMAL  DATA  FILES  FOR  PROM  3 
AND  PROM  4  PROGRAMMING 


02010101010101010101010101010101 
0101010101010.1010101010101010101 
01010101010101010101010101010101 
01010101010101010101010101010101 
01010101010101010101010101010101 
01010101010101010101010101010101 
0101 0101^101010101 01010101010101 
01 0101 010101 01010101 0.L01C  1010 101 
0404030303030303030303030303030.3 
03030303030303030303030303030303 
03030303030303030303030303020202 
02020202020202020202020202020202 
02020202020202020202020202020202 
02020202020202020202020202020202 
02020202020202020202020202020202 
02020202020202020202020202020202 
08080808070707070707070707070707 
070707070707060606060A0606060606 
06060606060606Q60606060606Q60605 
0 5 0 5 0 5 0 5 0 5 0 5 0 5 0 5  0  '5 0 5 0 5 0 5 0 5 0 5 0 5 0 5 
05050505050505050505050505050505 
05040404040404040404040404040404 
04040404040404040404040404040404 
04040404040404040404040404040404 
111010101010101.  OOFOF Of  OF OFOFOFOF  ' 
OFOEOEOEOEOEOEOEOEOEODODODODODOU 
ODODODODODODOCOCOCOCOCOCOCOCOCOC 
OCOCOCOBOBOBOBOBOBOBOBOBOBOBOBOB 
OBOBOAOAOAOAOAOAOAOAOAOAOAOAOAOA 
0A0A0A0A0A09 09090909090909090909 
09090909090909090909090908080808 
08080808080808080808080008080808 


BYTE 3  HEX 
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24242323232222222221212120202020 

1F1F1F1F1E1E1E1E1E1D1D1D1D 1C1C1C 

1C1C1B1B1B1B1B1B1A1A1A1.A1A1A1919 

19191919181818181818181717171717  4+ 

1 71 716161616161 61616 151 51 51515.15 

15151514141414141414141414131313 

13131313131313131212121212121212 

12121212111111111111111111111111 

5554535251 5050 4F4E4D4D4C4B4A4A49 

48484747464545444443424241414040 

3F3F3E3E3D3D3C3C3C3B3B3A3A393939 

38383737373636363535353434343333  ,+ 

33323232313131303030302F2F2F2E2E 

2E2E2D2D2D2D2C2C2C2C2B2B2B2B2B2A 

2A2A2A29292929292828282828272727 

27272626262626262525252525252424_ 

FFFBF7F3F0ECE9E5E2DFDCD9D6D4D1CE 

CCC9C7C4C2C0BDBBB9B7B5B3B1  AFAEiAC 

AAA8A6A5A3A 1 A09E9D9B9A9897959493 

91908FBE8C8B8A898887868483828180  , 

7F7E7D7C7B7A7A797877767574737372  L 

7170706F6E6D6D6C6B6A6A6968686766 

66656464636362616160605F5F5E5E5D 

5C5C5B5B5A5A59595B58575757565A55 

FFFFOOAA009955490038CCA2AAD82488 

0087 1CBC6618D1 90551 EECBD92694421 

00E0C3A68E755E48331F0CFAE8D8C8B9 

AA9C8F82766A5E53493E342B22191008 

OOF8FOE9E1DAD4CDC7COBAB4AFA9A49E  1 

99948F8A868 17D 7874 706 C6864605C58 

55514E4A4744413E3B3835322F202927 

2421 1F1C1 A1815131 10E0C0A08060402 


BYTE 3  HEX 
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FEFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFfrFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF  8 

FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FCFCFDFDFDFDFDFDFDFDFDFDFDFDFDFD' 
FDFDFDFDFDFDFDFDFDFDFDFDFDFDFEiFD 
FDFDFDFDFDFDFDFDFBFDFDFDFDFEFEFE 
FEFEFEFEFEFEFEFEFEFEFEFEFEFEFEFE  7 

FEFEFEFEFEFEFEFE F  E F E F E F E F E F  E F E F E 
FEFEFEFEFEFEFE F  E  F E F E F E F  EFEFEFEF E 
FEFEFEFEFEFEFEFE FEFEFE F E F  E F  E F E F  E 
FEFEFEFEFE  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E_ 

F  8  F  8  F  8  F  8  F  9  F  9  F  9  F  9  F  9  F  9  F  9  F  9  F  9  F  9  F  9  F  9  ‘ 
F9F9F9F9F9F9FA F A  F A  F A F  A F A  F  A F A F A  F  A 
F  A  F  A  F  A  F  AF  A  F  A  F  A  F  A  F  A  F  A  F  A  F  A  F  A  F  A  F  A  F  B 
FBFBFBFBFBFBFBFBFBFBFBFBFBFBFRFB  6- 

FBFBFBFBFBFBFBFBF8FBFBFBFBFBFBFB 
FBFCFCFCFCFC F C F C F C F C F  C F  C FCFCFCFC 
FCFCFCFCFCFCFCFCFCFCFCFCFCFCFCFC 
FCFCFCFCFCFCFCFCFCFCFCFCFCFCFCFC 
FOFOFOFOF  OF  OFOFOF IF 1 F 1 F 1 F 1 f 1F1F1 ' 

F 1 F2F2F2F2F2F2F2F2F2F3F3F3F3F3F3 

F3F3F3F3F3F3F4F4F4F4F4F4F4F4F4F4 

F4F4F4F5FSF5FSF5F5FSF5F5F5F5F5F5 

F5F5F<SF6F6F6F6F6F6FAF6F6F  6F6F6F6  5  ' 

F6F6F6F6F6F7F7F7F7F7F7F7F7F7F7F7 

F7F7F7F7F7F7F7F7F7F  7F7F7F  8F8F8F8 

F8F8F8F8F8F8F8F8F8F8F8F8F8F8F8F8 


BVT!i3  I ( I . X 
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DCDCDDDDDDDEDEDEDFDFDFBFEOEOEOEO 

E1E1E1E1E2E2E2E2E3E3E3E3E3E4E4E4 

E4E4E5E5E5E5E5E5E6E6E6E6E6E7E7E7 

E7E7E7E7E8E8E8E8E8E8E8E9E9E9E9E9 

E9E9EAEAEAEAEAEAEAEAEBEBEBEBEBEB 

EBEBEBECECECECECECECECECECEDEDED 

EDEDEDEDE D E D E D E D E E E E E  E E E E E E E E E  E E 

EEEEEEEFEFEFEFEFEFEFEFEFEFEFEFEF__ 

ACADADAEAFB0B1B1B2B3B4B4B5B6B6B7 

B8B8B9BABABBBBBCBDBDBEBEBFBFC0C0 

C1C1C2C2C3C3C4C4CSC5C5C6C6C7C7C7 

C8C8C9C9C9CACACACBCBCCCCCCCBCBCD 

CDCECECECFCFCFBODODOBOD 1 D 1  D:L  D2D2 

D2D2D3D3B3D3D4D4D4D4D5D5D5H5D6B6 

D6D6D&D7D7D7D7D7D8DSD8D8D8D9D9D9 

D9D9DADADADADADBDBDBDBDBDBDCBCDC 

03070B0F 1216191C :l F2225282B2E3033  * 

36383B3D3F414446484A4C4E50525455 

57595B5C5E5F6163646667686A6B6D6E 

6F 70727374 75777S797A7B7C7D7E7F 80 

81828384858687S8898A8A8B8C8D8E8F 

8F909 1929293949595969 7979899999 A 

9B9B9C9C9D9E9E9F9FA0A1 A1 A2A2A3A3 

A4A4A5A5A6A6A7A7A8A8A9A9AAAAABAB_ 

0000560067ABB700C8345E5628DC7800 

79E4449AE82F70ABE214436E97BCDF00 

203D58728BA2B8CDE1F406 1828384756 

64717E8A96A2ADB7C2CCD5DEE7F0F800 

0810171F262C333940464C51575C6267 

6C71767A7F8388809094989CA0A4A8AB 

AFB2B6B9BCBFC2C5C8CBCED1D4D7D9DC 

BFE1E4E6E8EBEDEFF2F4F6F8FAFCFE00 


BYTES  IlliX 
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00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000  8+ 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000  _ 
00000000000000000000000000000000  ' 

00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
000000000000000000000000000007KK7 
00000000000000000000000000000000 
00000000000000000000000000000000 
oooooooooooooooooooooooooooooooo  e,  + 

00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
OOOOGOOOCOOOCOOOOOOOOOOOOOOOOOOO 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 


BYTH4  III: X 
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00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000  4+ 

00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000  ,+ 

00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000  __  _ 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000  7  + 

00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
00000000000000000000000000000000 
7F7F  402A20191512 100E0C0E<0A090900 
0807070606 0605050505040404 040404 
04030303030303030303030202020202 
02020202020202020202020202020202  1  + 

020101010101010101010.1.010.10  10101 
0101010101010101010101010101010:1. 
0101010101010.1010101010101010101 
0101010101010:1.0101010101 01 0:1.0. i  01 


BYTI-4  HFX 
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FFFFFFFFFFF FFFFFFFFFFFFI- .  CFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFf  ITFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFF  Ff  FFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFF  FFFFF 
FFFFFFFFFFFFFFFFFFFFF FFF FFFFF FFF 
F F F F F F  F FF  F  F F  F  F FF  F  F F F  F  F  F I  ■  I  •  F  F F  F  F  F  F 
FFFFFFFFFFFFFFFFFF  FFFFF FFFFF FFFF 
FFFFFFFFFF F F F F F F  F F F  F  F F F F F F F F  FFF  F_ 
FFFFFFFFFF F F F F  F F  F F  F F F  F F  F  FFF  F  F F  F  F 
FFFFFFFFFFFFFFFFFFFFFFFFFFF FFFFF 
FFFFFFFFFF  Fr  FFF"  F  FFF  F  F  F  FFF  F  F  F  F  F  F  F  F 
FFFFF  FFFFFFF  F  F  F  F  F  F  F  F  F  F  F  F"  F  F  F  F  F  F  F  F 
FFFFFFFFFF FFFF F  F F F F F F  F  F F F F  F F F F  F F 
FFFFFFFFFFFFFFF F  F  F F  F  F  F  F F F  F F  F  F  F F F 
FFFFFFFFFFF  F F  F  F  F  F  F  F  F  F F F  F  F  F F F F F F F 
FFFFFFFFFFFFFFF  F F  F F F F F F  F  F  F F  F F F F  F 
FFFFFFFF FFFFFFFFFFFFFFF F F F F F F F F F 
FFFFFFFFFFFFFFFFFFF F F F  F  F  F  F F  F F F F  F 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFF FFFFF 
FFFFFFFFFFF  F F F  F F F  F F  F F  F F F  F F F F F F F F 
FFFFFFFFFFFFFFFFFFF7  F  F  F  F  F  F  F  F  F  F  F  F  F 
FFFFFFFFFFFFFFF FF F  F  F F F F F F  F F FFF  F F 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF" 
FFFFFFFFFFFFFFFFFFFFF F F F F F F F F F F F 
FFFFFFFF F  F  F F  F  F  F F F FFFFFFFF F F F F F F F 
FFFFFFFFFFFFFFFFFFFFFFFFFFF FFFFF 
FFFFFFFFFFF FFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFF FFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFF F  F F  FFF  F F F F  F F F F F F F 
FFFFFFFFFFFFFFFFFFFFFFFFFFF FFFFF 


BYTH4  MMX 
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FFFFFF F FFFF F F  F F F F F  F F F F F  F  F  F  F F F F F F 
FFFFFFFFFFFF F F F F  F F F F F F F F F F F F F F F  F 
FFFFFFFFFFFFFFFFFFFFFFF F F F F F F F F F 


FFFFFFFFFFFF  FFFFFFF FFFFFFf  FFFFF  F 
FFFFFFFFFFFFFFFFFF FFFFFFFFFFFFFF 
FFFFFFFFFFFFFF FFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFF'FFFFF  FFFF  FFFFFF  FFFFF 
FFFFFFF  F  F  F  F  F  F  F  F  F  F  F  F  F  F  FFFFFFF  FFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFF FFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFF FFFFFF 
FFFFFFFFF  FFFFFF  FFFFFFFFFFFF'F  FF  FF 
FFFFFFFFFFFFFFFFFFFFFFFFFF FFFFFF 
FFFFFFFFFFFFFF FFFFFFFFFFFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFF FFFFFFFFF 
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
FFFFFFF'FF'FF  FFFFF  FFFFFFF'FFFFFT  FFF 
FFFFFFFFFFFFFFFFFF FFFFFFFFFFFFFF 
FFFF'  F  F  F  F  F  F  f-  F  F  F  F  F  F  F  i::  F  F  F  F  F  F  F  F  F  F  F  F  F 
FFFFFFFFFFFFFFFFFF F'FFFFF'F"  FFFFFFF 
FFFFFFFFFFFFFFFFFF FFFFF FFFFFFFFF 
FFFFFFFFFFFFFFFFFF  FFFFFFFFFI-  FFFF 
FFFFFFFFFFFFFFFFFFFFFFF!"  FFI::  FFFFF 
FFFFFFFFFFFFFFFFFFFFFFF FFFFFFFFF 
80C0D5E0E6EAEDF0F1F3F4F5FAF6F7F8 
F8F8F9F9F9FAFAFAFAFBFBFBFBFBFBFC 
FCFCFCFCFCFCFCFCFCFCFDFDFDFDFDFB 
FDFDFDFDFBFDFDFDFDFDFDFBFDFDFDFE 
FEFEFEFEFEFEFEFEFEFEFEFEFEFEFEFE 
F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E  F  E 
FEFEFEFEFEFEFEFEFEFEFEFEFEFEFEFE 
FEFEFEFEFEFEFEFEFEFEFEFEFEFEFEFF 
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APPENDIX  F. 

FORTRAN  PROGRAM  FOR  CALCULATING 
EXTENDED-PRECISION  TABLE  ENTRIES 


INTEGER*2  YINT<2048) 

LOGICAL*!  YF'RGM < 32 ) 

EQUIVALENCE  <YINT( 1 ) »YPRQM< 1 ) ) 

REAL*16  Y ( 2048  ) fCF»C<16) »XI»XK»DEL»R 
C< 1 >=16128. 

C<2>=-192. 

C ( 3 ) =-4256 ♦ 

C<4>=-4240. 

C<5>=~3208. 

C ( 6 ) =-2180 . 

C ( 7 ) =-1410 • 

C<8)=~897. 

C( 9) =-114943 « 

C( 10) =-65727. 

C( 11 >=-37023. 

C( 12 >=-20623. 

C(13)=-11399. 

C( 14 >=-6275. 

C< 15 >=-3457. 

C<16>=-1920. 

CF=2  ♦ 

M=1 

DO  20  J=1 >16 

IF ( J  .GE.  9)  M=2 

L2=128*J 

L1=L2-127 

L3=8*M-J 

DO  20  I=L1>L2 

XI  =  I 

DEL=2.**L3 

XK=C<J>+XI*DEL 

IF( I .EQ.897  .OR.  I.EQ,898)G0  TO  15 
IF< I .GT .897  .AND.  I.LE.1024>G0  TO  10 
IF( I .GE. 1921  .AND.  I.LE.2048>G0  TO  40 
IF ( I .EQ . 128  .OR.  I. EQ. 11 52) GO  TO  30 
IF (M.EQ«2)CF=~2« 

Y(I)=CF/(2.*XK+DEL-1. > 

GO  TO  20 
10  Y(I>=1./XK 
15  Y<897>=1.~(2.**(-15> > 

Y<898)=Y<897) 

GO  TO  20 


TABLE2  FORTRAN 


40  Y(I)=-1./XK 
GO  TO  20 

30  Y(I)=CF/(2.*XK+255. ) 

20  CONTINUE 

DO  70  1=1 » 2048 
R=Y  < I ) 

N=1 

IF (R.GE.O)GO  TO  45 

N=-l 

R=~Y < I ) 

45  R=R-2.**(~15) 

IF'  ( R ,  LT ♦  0  > GO  TO  50 
GO  TO  45 

50  Y(I)=R+2.**<  -15) 

IF ( N . EQ ♦ -  1 ) Y  < I ) =-Y ( 1 ) 

WRITE(3»3)Y(I) 

3  FORMAT (F36. 33) 

N=1 

ICUM=0 

R=Y(I) 

IF  <  Y ( I ) . GE . 0 ) GO  TO  25 
N— 1 
R=-Y ( I ) 

25  R=R-2.**<-30) 

IF ( R . LT .0) GO  TO  35 
ICUM=ICUM+1 
GO  TO  25 

35  YINT(I )=N*ICUM 
70  UR I TE  <  4 » 4 ) Y I NT  < I ) 

4  FORMAT (18) 

DO  60  1=1 f 4065*32 

WRITE ( 1 » 2 )  YPROM ( I ) *  YPROM ( I  +  2 ) , YPROM ( 1+4 ) > YPROM ( 1+6 >  r 
1YF'R0M<  1+8) » YPRQM<  1  +  10  ) » YPROMC 1  +  12)  >YF'ROM(  1  +  14)  r 
1 YPROM ( I  + 16  > » YPROM  ( 1+18 ) » YPROM  < I +20 ) » YPROM ( I +22  >  r 
1 YPROM ( I +24 ) , YPROM  ( I +26  >  *  YPROM  < I +28 ) r  YPROM ( I +30 ) 

60  WRITE (2» 2)  YPROM  <  1  +  1 )  fYF'ROM<  1+3 )  r  YPROM  ( 1+5)  t  YPROM(  1+7 )  > 
1 YPROM ( I +9 ) r YPROM < I + 1 1 ) f YPROM ( I + 1 3 ) > YPROM ( I + 1 5 ) » 

1  YPROM  ( 1+1 7 ) > YPROM  <  1  +  1 9 )  »  YPROM  ( 1+21 ) » YPROM  ( 1+23 )  r 
1  YF'ROM  (1+25)  r  YPROM  ( 1  +  27 )  t  YPROM  ( 1+29)  >  YF'ROM<  1  +  31 ) 

2  FORMAT (32Z2) 

STOP 

END 


TABLFZ  FORTRAN 
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APPENDIX  F 

HEXADECIMAL  DATA  FILES  FOR  PROM  1 
AND  PROM  2  PROGRAMMING 


03020A1A315075A2DA105199E73B95F5 

5AC637AD29AA30BB4BE07AL8BB620EBE 

722AE7A76B33FFCFA2795331 12F6DEC9 

B7A89C948E8B8B8D939BA6B4C4D6EB03  8+ 

1D3A58799DC2EA14406F9FD2063C75AF 

EB2969ABEF347BC40E5BA9F8499CF046 

9DF650AC0968C8298CF056BC258EF964 

n240AF2092057AEF65DBS6CF4AC64381 

1 AO A09 193866 A3EF48B025A735D 1792B 

EDB88F715E5557647A9AC4F73379C81F 

7FE758D051D96901A047F4A96426EFBF 

9470533B2A1E 18181 D28394F6A8AAFD9  7+ 

083D75B3F53C87D72B83E040A50E7BEC 

60B955D558DF6AF88A  .1.  F  B 7 5 3 F 2 9 4 3 9 E  2 

8D3CEEA25AJ.4B191.541AE2AD7B4Bi.BF3 

CBA5826142260CF5DFCCBCADA1968E88 _ 

EB874524234280DD56EDA06E585B79B0 

0067E77E2BEFC8B7BBD3004093F972FD 

9A4908D9BAABACBDDD0C4B97F25BD2S6 

E88632EAAE7F5C45393843597AA5DB1B  6+ 

66BA1981F26DF27F16B55D0DC6885123 

FCDEC7B8B0B0B7C5DAF71A4475ACEA2E 

79C9217EE14AB92EA828AE39C95FFA9B 

40EB9A4F08C68951 1DEEC49E7C5F4631 _ 

081E7F291A51CB8987C43FF6E9157A16 
E8EF29963504022F8910C29FA6D62EAE 
54201 12760BC3BDB9D7F81 A2E240BB54 
0ADBC8D0F33086F67F20D9AA9291A6D2  5+ 

•  1 369D555E A924F1F02F7001 A4785D535 
A72ABCSF12D5A788797886A2CC054B9F 
006EEA7207A95711D7A9877065657086 
A7D2084993E847B0239F2SB44CEE994D 


BYTE!  HEX 
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1 


EDD2AFBFBDAAB4A43259]  IAA543AC9E  1  4 
0C8275E0C217DD11B2BC2E053FDAB42A 
DCE748000BA81 A135DF  2D2FBAB22 1 DSC 
DD9FA1 E15F 1 90F3EA74921 2F73EB9775  4 + 

85CA37D7AAA2CB20A14D23224A9A1  2B0 
755FAFA3FB7614D5B7BBE0258A0EB173 
5350ABA2FAAAF2985A3A2C3CA5A8037A 
02A5A032 1 A 1 A2F5 A9BF25DDD72 j  CB9AA_ 
2A6EF7B5958980AB3CEA5A8C71FB1FH2 
OBBDEOAA51 8E 1 AE 1 E9248C 1 8C3845A32 
11EFC48C41DC5BB6EAF2C86AD1FBE38S 
BDE7A0041 1C1 12028CAEA5AD85E9D74C  -  + 

46C2BE372B987CD49FBA849A1C055AOC 
2AA27DB84F428F3431. 8329 2 2 A  D  0  8  F  3  2  B 
AF809AFBA99BD3501 1 1459DFA5AAECAC 
281F51BCA13D51 9B1 ACFB8D524AAS93D 
7  F  4  2 1 F  A  A  8  0  C  8  B  2  B  A 5  3 0 F  7 8 2 1 A 1 9  9  A  A  7  C 
BC194903030ACE20C380247C5A90F358 
978B0EFE389C0BAA906FEADD3AEAC9CD 
BEEAD08B0428EA2EEE199EAE7CBA1A90  7+ 
0F8CF94B7B793EBEF0CB455AF41 8B8CD 
5039801E0C45C0776583CC39C5AB25ED 
C098AF420CC8720A7FDA12240CCA4EA1 

BC9A3A97AF7E013A19A8E0BE3F6223EQ _ 

0000005500CCAA4900C7AASD55C42444 
00C3E3D733C32E0BAA28A297922C2221 
00 1 FE 1 4 1 7 1 9 .1 AB4 1 99 18A1 F  4 .1 7  Cl  8582 
550A1 44 1 B 1 CF4B1 2494796E5 1 1 .1 4 1041 
00C00F8DF003A0B43838C80D357B2074  1 + 

CC870CC830C0FAA48B02A040420B4.1. 91 
AA3F05B50AC3A065D8C0E71A25DA0985 
24BC2335CBC0F24O88AB8A0808AE2004 


HYTlil  HEX 


55 


0201FAEAD3B48E612EF3B2AB1BC96F0F 
A93DCC5ADA59D348B82389EB48A1F545 
91D81C5C97CF0334A08AAFD2FOOC2439 
4B5A666F747777756F675C4E3E2C16FF 
E5C8A988653F17EDC1936230FCC58D52 
1 6D89856 1 3CD863DF3A75909B86 A 1 1 BC 
A40BB155F89939D8751.  1 AC45DD73099D 
30C152E1 AFFC88 129C24AB32B73BBE81. 
F60706F7D7A9AB20CA5FEAA7D83C94E0 
20547D9BAEB7B5A892724814B89243EC 
8C23B33AB930A008A9C315A0A5E31 A4A 
7498BACDDFEAF0F0EADFCF  B99E7E582E 
FFCB925412CB8030DC8327CAA1F88B1B 
AA2DB131 AE279C0E7CE74EB31471CC24 
78C918A3ABF 13474B1EB23588ABAE712 
3AA083A4C2DE.F81025384957A4AE7A7C 
59BCFD1D1CFCBD60E54E9ACBE1DCBE8A 
35CD4DB50743A879755C2FEE9A33BA2E 
91E221506F  7D7CAA4A1ADB8E33C952CE 
3C9DF13873A2C5BCE7E8DCCAA5794302 
B8A3049C2AAF2A9C0AAABE0D5492C8F  7 
1D3B52A1A8A8A1523D20FDD3A2AA2CE8 
9D4CF59734CB5CE7ACECAADB4AB41978 
D32878C40A4C89C1F  4234E7495B2t;BE0 
19FE99EAF5BA3B79783AB7FC0ADAAECE 
F9EFB141 A0CECC9D40BA012218EA8B09 
A1929F874BEDACCA0A221EFBB959DC41. 
8AB7C8BE9A5C039208A5AAD8EFEED8AB 
681 1A4228CE2255470796F542AE 79634 
C13EAA06528EBBD9E8E8B9BC9:L5710BC 
5 AEAAEE54F ABFE437CA9CBE 1 ECEBDFC9 
A77B4403B8A202992AA92393F A57ACF 8 


BYTIi  1  Ili:X 
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4044912B1759F7F3521848E7F77D7BF4 
EC6664EAFA97C380D2BB3D59136C6705 
4832C503EEB6CFC875ri6EDBB42827F38 
AEE4DB930I.I4B4F18A7FF200BC14290AB 
954ED8325D5B2DD24B9ABFBB8E39BD19 
506 1 4B1 4B83895D1 EAE2B970087F D8  1 3 
2F2E10D57D0A7BI»10C2I.»3320F4AE50DA 
4BA5E71227250DDF9B42D4S1BA0E4E7B 
FD260DC65FE9740EC5A8C323D4E25944 
AC9C1 F3DFFAF95 7A259EEC  1. 8272 10CEE 
CFB3A09DAED 9239128 E  C  E  2  0  F  7  A  .1.  D  0  A  3  A 
B 1 7B9708D3F A806AB A72972A2F A89802 
E74B309885FAF9849D46B04FB4B04A76 
44B0BDABBBB35095831C6053F34447F  B 
648154DE21 1DD44674610B75A08B39AA 
E0r*A99206D83620A7riBBC59C40B2F202 
8E27FFBBF94E47672C0B75D40CFC806D 
1 4  C  4  C 6 6 2 B  8  6  B  5  7  D  5 1 F  6  8  E  4  C 2 3 3 6 1 7 9 A  2 
03C305EC992CC37D75C78DE0D88E188B 
FC802B0F3ECBC64049EF43522AD96DF0 
71FB99573F5CB95E57ABA58C2A47EB1E 
E84F5C1 480 A6SC38B.IFD  2224  0E<BA984A 
F59C46F8B480615A6FA500833 110226 A 
EDAEAFF480577AEEB5D 1461 A 4 3D 1 C218 
0000AB0O3456B7OO399AA3AB3CDCBCO0 
3  D 1 D  2  9  C  D  3  D  D  2  F 5 5 6 D  8  9  E  6  9  A  E  D  4  D  E  D  F  0  0 
E 1 1FBF8F  AF"95BF  67E89F  0CE93F7B7EAB 
F6ECBF4F31B5EEB7B96A1 BEFECFOBFOO 
40F17310FD604CC8C838F3CB85E08C34 
79F438D040069C7SFEA0C0BEF5BF6F56 
C1FB4BFA3D609B284019E6DB2AF77BDC 
44DDCB35400EC0785576F8F892EOFCOO 
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017F7D7B7977757371706E6C6A696765 
6462615F5E5C5B59585655545251504E 
4D4C4A4948474544434241403F3D3C3B 
3A393837363534333231302F2E2B2C2C 
2B2A2928272625252423222121201F1E 
1D1D1C1B1A1A191 818171 61 515141313 
121111 1O1OOFOE0EODOC0COB0B0AO9O9 
08080707060605040403030202010100 

06027E7A76726E6A6763605C5955524F - 

4B4845423F3C393633302D2A28252220 
1D1A181513100E0C09070402007E7B79 
7775737 16F6D6B696 7 65636 15F5D5B59 
58565452504F4B4B4A4846454342403E 
3D3B3A383735343231302E2D2B2A2927 
2625232221201E1D1C1B191817161513 
121 1 100FOEODOCOA0908070605040302 
1C140C047C746C645D554E474039322B 
251E171 10B047E78726C67615B55504A 
45403B35302B2621 1C18130E0905007C 
77736F6A66625E5A56524E4A46423E3B 
3733302C2825211E1B1714110D0A0704 
007D7A7774716E6B686563605D5A5755 
524F4D4A474542403D3B383633312E2C 
2A272523211E1C1A1815131 10F0D0B09 
006E5C4B3A291808786859493A2C1D0F 
007265574 A3B3023160A7D7165594E42 
372C21 160B00766B6 157 4D 43393026 ID 
140A01786F675E554D453C342C241C14 
0P057D766E676059524A443D362F2822 
1B150E08027B756F69635D57514C4640 
3B352F2A251F1A150FOA05007B76716C 
67625E59544F4B46423D3934302B2723 


BYTE 2  MiX 
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360D6531H6704A25015D3A 17755332 12 
7252331 476593 B1F02664B301 57 A6047 
2D147C644C341D066F58422C17026D58 
432F1B0774614E3B28160472604E3D2C 
1B0A7AA95949392A1A0B7C6D5E4F4132 
2416087A6C5F5244372A1D1 104786B5F 
53473B2F23180C01 766 B60554A3F352A 
20 150B01776B63594F453C 322920 160D 
001022374E67021F3E5E01.254B721C4Z" 
73204F0032651 A4F063F78336E2B6928 
68286A2D71357B41085019632D784411 
5E2C7B4B1B6B3BOF6134085C31065C33 
0A6139126B441E78532E0A6643207D5B 
39 17765535 1575563 7 187A5C3E2 10367 
4  A2E12765B40250B71 5731*230 A7 15840 
28 10 78604932 1B046E57412B16006B56 
000C2750074B1C796256545D6F0B305E 
145319673C187A64534945464D596A01 
1 C3C61 093768 1E57 14551 96 12D7B4D22 
7A553314785E46321F1002776E676360 
6061646A717A051 1 1F2F4154687F162F 
4A660322426305294E741B446D184470 
1E4D7D2E6012467B30661E560F48033E 
7A3775337232733476387B3F04490F55 
00000055004C2A1200716674554E4944 
0043385033303A592A5C275024772204 
007C61 1D1C4F281A191C18055D162C18 
555E2E411321686412707B31 11 174210 
000F3E0970730E420E7067741 44701*64 
4C454E650C4002522E170B0B162C4C76 
2A682F7E5737201 10909101 F344F721 A 
497D38783D08582E08674B3421 120802 


BYTI-2  1 1 1: X 
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■  V- 


7 


FF81828486888A8C8E8F91939596989A 
9B9D9EA0A1A3A4AAA7A9AAABADAEAFB1 
B2B3B5B6B7B8BABBBCBDBEBFC0C2C3C4 
C5C6C7C8C9CACBCCCtiCECFD0DiD2D3D3 
D4D5D&D7D8D9DADADBDCDDDEDEDFE0E1 
E2E2E3E4E5E5EAE7E7E8E9EAEAEBECEC 
EDEEEEEFEFFOF IF IF2F3F3F  4F4F5FAF6 
F7F7F8F8F9F9FAFBFBFCFCFDFDFEFEFF 
F9FE8285898D91 95989C9FA3AAAAAITF0 
B4B7BABDC0C3C6C9CCCFD2D5Li7DADBriF 
E2E5E7EAECEFF1F4FAF  8FBFDFF  81848A 
888A8C8E90929496989A9C9EA0A2A4AA 
A7A9ABADAFB0B2B4B5B7B9BABCBBBFCJ. 
C2C4C5C7C8CACBCDCECFB1D2D4D5DADB 
D9DADCDDDEDFE1E2E3E4E6E7E8E9EAEC 
EDEEEFFOF 1F2F3F5F6F7F  8F9FAF  BFCFD 
E3EBF3FC848B939BA2AAB1B8BFC&CDD4 
DBE1E8EEF5FB81878H93999EA4AAAFB5 
BABFC5CACFD4D9HEE3E8ECF1F  6FAFF  83 
888C9095999DA1  A5A9ABB.1B5B9BDC1C5 
C8CCB0D3D7DADEE1E5E8EBEFF2F5F8FB 
FF8285888B8E9194979A9C9FA2A5A8AA 
ADB0B2B5B8BABDBFC2C4C7C9CCCED1D3 
D5D8riAnCDFElE3E5E7EAECEEF  OF  2F4FA 
8 1 92A4B5C6D7E8F88898A7BAC6D4E3F 1 
FF8D9BA9B6C3D0DDEAF6838F9BAAB2BE 
C9D4DFEAF5FF8A949FA9B3BCC6D0D9E3 
ECF5FE879099A2AAB3BBC3CBD3BBE3EB 
F3FB828A9198A0A7AEB5BCC3CAD0D7DE 
E4EBF1F8FE848A909A9CA2A8AEB4BABF 
C5CAD0D5DBE0E5EBF0F5FAFF  84898E93 
989DA2A7ABB0B5B9BEC2C7CBCFD4D8DC 


BYTE 2  HEX 


C£F79FC7£E94B9DE83A7CAEC8EB0n  1 F 1 

91B1D0EE8CAAC7E4809CB8D3EE88A2BC 

D5EE869FBACEE5F  C93A9BFD5EB8095AA 

BED2E4FA8EA1B4C7D9EBFE90A1B3C4D5 

E6F78798A8B8C8D7E7F A8594A3B2C0CF 

DPEBF98794A2AF  BCC9DAE3F0FD8995A2 

AEBAC&D1DDE9F4FF8B96A1ACBAC1CCD6 

E1EBF5808A949EA7B1BBC4CED7E1EAOL 

9B8BF8E2CBB 1 9AF9D9B895F 1CAA2F9CE 

A1F3C493E0ADF8C28BD298DEA2E5A7E7 

A7E6A4E19DH893CC85BCF3AADF94C8FB 

ADDF90C1F09FCEFCA9DA82AED983ADD7 

FFA8D0F79EC4EA90B5DAFEA2C5E88BAD 

CFF091B2ri2F292BlD0EF8I.iABC8E6839F 

BCD8F48FABC4E0FB95AF  C9E2FB94ADC5 

DDF58DA5BCri3EA8J97ADC3It9EF8499AF_ 

FBE8C595D78DBAD3E5ECE8DAC3A2F9C7 

8DCA80AFDAF791  A4B2B9BABAAD9E8AF.1. 

D4B18BDFB0FriC58ACB88C2F8ABI.iB88BJ 

D7FB9CBAn5ED8397A8B6C3CDD5DADEDF 

DFDCD8D2CAC0B4A79887F5E1CCB59C83 

E7CBAri8EEDCBA884DEB790E7BH91E5B8 

89DAAAF8C693DFAAF4BB8ACD94DA9FE3 

A6E9ABECADEDACEAA8E5A2DE99D3SDC7 

OOOOAAOOB3D5ED008E998BAAB1B6BBOO 

BCC7AFCCCFC5AAriSA3D8AFDB88LiDFB00 

839EE2E3B0D7E5E6E3E7FAA2E9B3E7AA 

A1IUBEECDE979BEB8F84CEEEE8BDEF00 

F0C1F68F8CF1BDF.18F988BEBB8F29BB3 

BAB19AF3BFFDADDlE8F4F4E9D3B389ri5 

97rt081A8C8DFEEFAFAEFE0CB&08riE5BA 

82C787C2F7A7MF798B4CBDEEDF7FB00 


BYTE2  HEX 
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APPENDIX  G 

SAMPLE  CALCULATIONS  FOR  TWO  PSEUDODIVIDE  TEST  CASES 

Test  Case  1 


X  = 

3321 

(integer) 

=  CF9 

=  0.101 

348 

876  953  125 

Y  = 

1672 

(integer) 

=  688 

=  0.051 

025 

390  625 

The 

exact 

decimal 

value 

for  1672/3321 

=  Q  = 

0. S03  462  812  405  901  ... 

R  =  2'15/X  =  0.000  301  108  ... 

A.  A 

Q  =  R-Y  is  to  be  calculated  in  two  parts,  with 

R  =  R  +  R  ,  and  Q  =  Q  «■  Q  «  R  * Y  +  R, *Y  . 
m  1  m  xl  m  1 

A 

Some  scaling  maneuvers  are  needed.  If  we  interpret  the  Rm 

- 1 5 

table  as  being  scaled  by  2  ,  but  the  R^  table  as  not  scaled, 

A 

then  Qm  can  be  retrieved  from  the  multiplier  LSP  with  assigned 
- 1  - 1 5  ~ 

weights  of  2  to  2  ,  and  may  be  found  subsequently  at  the 

multiplier  MSP  with  weights  interpreted  as  2  ^  to  2 

Integer  X  =  3321  is  the  88th  location  in  reciprocal  zone  5+, 
where  quantization  factor  (A)  is  16  and  integer  X  ranges  from 

A  A 

1920  to  3967.  From  the  hex  tables,  R  =  0009  and  R.  =  6F88. 

m  1 

Accordingly , 

Qm  =  (0009)  (688)  =  3AC8  =  0.459  228  515  625 

Q:  =  (6F88) (688)  =  2D87040  =  05B0  =  0.044  433  593  750 


62 


and 


Q  =  Qm  +  Qi  =  0.503  622  109  375 
Error  ET  =  (Q  -  Q)/Q 

=  (0.503  662  109  375  -  0.503  462  812  405)/ 
0.503  462  812  405 
=  0.00039 


Test  Case  2 

X  =  -3321  (integer)  =  F307  =  -0.101  348  876  953  125 
Y  =  1672  (integer)  =  688  =  0.051  025  390  625 
The  exact  decimal  value  for  1672  /  (-  3321)  =  Q  = 

-0.503  462  812  405  901  ... 

R  =  2"15/X  =  -0.000  301  108  ... 

Integer  X  =  -3321  points  to  reciprocal  zone  5*,  88th  location. 

A  A 

From  the  hex  tables,  R  =  FFF7  and  R,  =  90D9. 

’  m  1 

Thus 

Qm  =  (FFF7) (688)  =  687C53S  C538  (selecting  LSP) 

=  -0.459  228  515  625 

Qx  =  (90D9) (688)  =  D2A0948=>  FA54  (selecting  MSP) 

=  -0.044  311  523  437  500 

and 

Q  =  Q  +  Q  =  -0.503  540  039  062  5 
1 

Error  ET  =  (Q  -  Q)/Q 

=  (-0.5035400390625  +  0.5034628124059)/ 
(-0.5034628124059) 

=  0.00015 
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UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (then  Data  Entered) 


()  7  Report  documentation  page 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

t  REPORT  NUAnnTft 

UW  KSI^Fr-^W 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT  S  CATALOG  NUMBER 

4-  TITLE  Submit*  > 

A  Compact  liartlware  Realization  for  Approximate  Division  . 

OF  REPORT  &  P  Eftfd&^C*E  RED 

^  (y  Technical  jfol J 

6.  PERFORMING  ORG.  REPpUA  iMIMCt-Rf 
Technical  Note2l979-57  j 

!  (W  .  i  f<r 

j  '  "yUbert  11  Jluntoon 

j  i - - - — - - J  t 

8.  CONTRACT  OR  GRANT  NUMBE  Rt, ) 

'  F19628-80-C-0002  f 

|  9.  PERFORMING  organization  name  and  address 

i  Lincoln  Laboratory,  M.l.T. 

:  P.O.  Box  73 

Lexington,  MA  02173 
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